JOURNAL OF RESEARCH of the National Bureau of Standards 
Vol. 87, No. 4, July-August 1982 

Measurements of the Octanol /Water Partition 
Coefficient by Chromatographic Methods 

Stanley P. Wasfk,* Yodu B. Tewari,* Michele M. Miller* 

National Bureau of Standards, Washington, DC 20234 

and 

J. H. Purnell 

Department of Chemistry, University College of Swansea, Swansea, Wales, U.K. 

June 3,1962 

A theoretical relationship is developed to provide a quantitative definition of hydrophobicity using 
established theoretical and semi-empirical relationships. A method of predicting partition coefficients of 
relatively water-insoluble third components between water and an immiscible second component is devised 
and tested. Comparison with experimental data for four classes of compounds in the water/n-octanol 
system at 25° C shows excellent agreement, indicating that values for substances for which direct deter- 
mination is experimentally precluded can be calculated with confidence. 
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1. Introduction 

In recent years there has been an increased interest in 
the use of hydrophobic parameters to study the fate of. 
toxic substances in the marine environment since the 
ability of organic compounds to bioconcentrate is be- 
lieved to depend upon the partition behaviour of 
molecules between lipid and aqueous phases [1,2]A An 
important and simplifying observation has been that by 
Neely, Branson and Blau [1] who demonstrated that 
bioconcentration factors for chlorobenzenes and 
chlorophenols between trout muscles and dilute solu- 
tions in water could be successfully correlated with then- 
partition coefficients in the n-octanol/water system, 
K o/w . Subsequently, Dunn and Hansch [3] compiled 
hydrophobic interaction data for a large number of 
organic compounds and showed that these could, indeed, 
be quantitatively correlated with partition coefficients of 
organic/ water systems. The weight of evidence has led 
Leo [4| to suggest not only that hydrophobicity is the 
most important parameter in bioaccumulation and 
biotransport but that this can be confidently determined 
in terms of octanol/ water partition coefficients. 

Because of the thousands of compounds being studied 
as potential hazards to the environment through bioac- 



cumulation, simple economics makes it desirable to 
devise some system whereby we measure values of K / w 
for key compounds which may be used to calculate 
values for related compounds. In the only approach to 
date, Hansch, Quinlan and Lawrence [5] have developed 
a method for estimating K / w based on additive group 
contributions; these group contributions or "11 values" 
being as defined by eq ( 1), 



n x = bgKj /w -logKfc 



o/w 



<1) 



♦Center for Chemical Physics, National Measurement Laboratory. 
'Figures in brackets indicated literature references at the end of this paper. 



where KJ/ W and K h /w are the octanol/ water partition 
coefficients for the derivative and the parent compound, 
respectively. Not surprisingly, in the light of the 
numerous such correlations established in GLC studies 
[6], n values rre often additive, and the method has met 
with acceptance. However, because of steric, electronic 
and hydrogen bonding effects there are many series of 
compounds for which the method fails. Some alternative 
would, therefore, be useful, particularly if the correlation 
method involved real and measurable physical proper- 
ties of the molecules concerned rather than a purely em- 
pirical set of parameters. This paper outlines one such 
approach and an indication of its applicability. 

Partition coefficients are generally determined by 
some variant on the traditional shake-flask method. This 
method is slow, tedious, often wasteful, and demanding 
in the standard of purity of materials it requires. Conse- 
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quently there have been many attempts to develop 
chromatographic methods to which, in principle, and 
normally in practice, none of the above objections apply. 
In the method presented in this paper K 0/ , w is defined, as 
proved later, by eq (2), 



Kq/w — 



= yff$ 



(2) 



where yY and y? are the activity coefficients at infinite 
dilution, based on volume fraction, for the solute in 
water and n-octanol, respectively. The quantity, yY t is 
calculated from the solute aqueous solubility, C w , and 
the solute molar volume, V\l\, and y? is determined from 
the corrected retention volume of the solute eluting from 
a column containing n-octanol as the stationary phase. 
Thus K o/w (at infinite dilution) can be calculated from yV 
and y?, which can be measured by two independent 
methods having all the advantages of the 
chromatographic approach, while the recognition of eq 
(2) opens the route to the alternative approach to be 
described. 

2. Theoretical 

The octanol/ water partition coefficient, K o/w , is 
defined as the equilibrium ratio of the molar concentra- 
tion of solute x in octanol, C , and the concentration in 
water, C w , in an octanol/ water system, viz. 



K o/w — 



(3) 



But, self evidently, K o/w may also be defined by eq (4) 

o/w 



K 



(4) 



a/v 



where K a / is the air/octanol (saturated with water) and 
K a/w is the air/water (saturated with octanol) partition 
coefficient. GLC theory yields the expression for K a / Q 



K a/o — 



= CRT 



nJtT 



= (V X +V +V W )P 



(5) 



where n x is the number of moles of solute x in the octanol 
(saturated with water), V x , V and V w are the volumes of 
solute, octanol and water, in the octanol phase, respec- 
tively, R is the gas constant, and P is the partial pressure 
of solute x above the solution at temperature, T. 

For an ideal vapour the ratio P/P° t where P° is the 
solute saturation vapour pressure, is equal to the activity 
(a). We may define the activity in terms of any quantity 
that defines relative amount, and an appropriate activity 



coefficient, e.g., we may write, 

a =r I x= y+O = y w w - y c C/C e 



(6) 



where x represents mole fraction, ^ the volume fraction, 
W the weight fraction, and C/C° the ratio of concentra- 
tion in solution to that in the pure solute liquid. 
Although the first of these definitions is the one most 
widely used we choose the second for reasons that will 
emerge. 
Thus, setting 



P 

substitution for P in eq (5) yields 
■K„/„ = 



(7) 



/a 



RI _ RT 



v r +v n +v» 



yfypo V(x)yIP° (8) 



where y? is the activity coefficient in water saturated 
with n-octanoL Provided there is no significant excess 
volume of mixing in either solvent, the normal solution, 



^a/w — 



RT 

VMy^Po 



(9) 



where yY is the activity coefficient in water saturaterd 
with n-octanol. Provided there is no significant excess 
volume of mixing in either solvent, the normal solution, 



K o/w — 



xy/y? 



(10) 



It remains only to emphasize again that y? and yY are 
values. An immediate and obvious attraction of eq (10) is 
the absence of explicit solvent parameters, which is not 
the case if y x is used. 

3. Dependence of the thermodynamic func- 
tions associated with solute partitioning be- 
tween an organic and an aqueous phase on 
vapour pressure. 



Hoare and Pumell [8] have shown that for solutes of 
similar chemical structure the GLC specific retention 
volume of solute x, VJx), is related to the saturation 
vapour pressure of the solute. P°, by the expression 



log V°{x) = -a log P° + constant 



(11) 



where a is a series constant. The validity of eq (11) has 
subsequently been further established for a wide range of 
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chemical types of both solvent and solute [9, 10, 1 1] and 
in a summary by Purnell [8]. In bo far as data are 
available there seems to be no recorded exception to the 
rule. Further since 

V°(x) = ^ 
e y + P° V(x)q 

where q is the solvent density, eq (1 1) yields 

log y? V{x) = (a° - 1 ) log P° + constant (12) 

where a° is the slope of the log V°{x) vs log P° plot for 
each series of solutes in n-octanol. 
In similar manner we obtain 

logyrHx) = (a w -l)logP° + constant (13) 

where a w is the slope of the log V (x) vs log P° plot for 
each series of solutes in water. Combining eqs (10), (12), 
and (13) we obtain 

log K o/vf = (a - a w ) log P° + constant (14) 

Thus for solutes of similar chemical structure the ther- 
modynamic functions, log P°, log y?K(x) and log yj^fx), 
associated with the partitioning of a solute between an 
organic and aqueous phase may each be expressed in the 
form 



Table 1. Solute Activity Coefficients and Octanol/ Water Partition 
Coefficients at 25.0°C. 



function = B + A log P 



(15) 



where B and A are numerically defined by the function 
being considered. 

4. Results 

The solute activity coefficients in octanol, yg and in 
water yV listed in columns 1 and 2 of table 1 were ob- 
tained from an earlier publication of Wasik et al (12). yj 
were calculated from the solute specific retention volume 
Vg, obtained from retention times of solutes eluting from 
a GC column containing n-octanol as the stationary 

phase, y^ values were calculated from solubility data us- 
ing the following equation 



f= l 

V 

iw 



(16) 



where <|> w is the solute aqueous solubility in volume frac- 
tion. The last two columns show a comparison between 
the calculated K / w (using eq (10)) and the experimental 
values. The agreement between the two sets of data is ex- 
cellent. 





LogrJ 


LogrJ' 


LogK o/w 


Solute 


Calculated" 


Experimental 




{Generator 








(GLC) 


column 
/HPLGorGC) 






n-Pentane 


0.555 


4.19 


3.63 (3.68) 


3.62 


n-Hexane 


.530 


4.73 


4.20(4.22) 


4.11 


□-Heptane 


.517 


5.28 


4.76 (4.77) 


4.66 


n-Octane 


.512 


5.80 


5.29 (5.29) 


5.18 


1-Hexene 


.504 


3.98 


3.48 (3.47) 


3.39 


1-Heptene 


.491 


4.58 


4.09 (3.07) 


3.99 


1-Octene 


.479 


5.24 


4.76 (4.73) 


4.57 


1-Nonene 


.470 


5.81 


5.34(5.30) 


5.15 


Toluene 


.509 


3.17 


2.66 (2.66) 


2.65 


Ethylbenzene 


.505 


3.66 


3.15(3.15) 


3.13 


n-Propylbenzene 


.494 


4.22 


3.73 (3.71) 


3.69 


Ethylacetate 


.621 


1.15 


0.53 (0.64) 


0.68 


n-Propylacetate 


.534 


1.64 


1.11(1.13) 


1.24 


n-Butylacetate 


.425 


2.15 


1.73 (1.64) 


1.82 



'Values in parentheses were calculated using the hydrophobicity 
equation (eq 17). 

5. Discussion 

Our approach of defining K o/w in terms of volume 
fraction based activity coefficients has several important 
consequences. First, the infinite dilution coefficients log 
K o/w , log y?, and log y£ for a given solute type are clear- 
ly described by a Knear relationship with the solute 
saturation vapour pressure. 

The extent to which this is true may be gauged by con- 
sideration of the values of the correlation coefficients, r, 
listed in table 2 derived from data of table 1 via linear 
regression as a consequence of the above. 

The data in table 1 indicate that values of log y% are 
large and change rapidly with log P°; whereas, log y , 
are much smaller and remain fairly constant (0.510 ± 
0.045) for all compounds. Thus the relevant values of log 
K o/w are determined by log yV. Our approach leads us to 
a quantitative theoretical definition of the hitherto em- 
pirical concept of hydrophobicity, H, 

H = \ogK o/yf = log y£- 0.510 

= logyy-fc 

The numerical constant k will vary for different solvents 
which means that H can be defined with respect to 
solvents other than n-octanol. The values of 
hydrophobicity {Hi calculated using eq (17) are listed in 
table 1 in parenthesis. A quantitative definition of H is 
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Table 2. Coefficients of the Regression Equation and the Coefficient 
of Correlation for the Solute 



Type of 
Correlation 


Type of 

Solute 


Slope 


Inter- 
cept 


Coefficient of 
Correlation 

(r) 


AbgP + B 












Alkanes 
Alkenes 
Aromatics 
Acetates 


-0.9957 
-1.1400 
-0.8809 
-0.7700 


6.298 
5.981 
3.770 
2.484 


0.999 
.999 
.998 
.999 


fe l n c+ fe 2 












Alkanes 
Alkenes 
Aromatics 
Acetates 


0.5230 
.5860 
.5200 

.5700 


0.993 
-0.1200 

2.117 
-0.4633 


0.999 
.999 
.999 
.999 


logK o/w = 

w+ki 












Alkanes 
Alkenes 
Aromatics 
Acetates 


0.01748 
.02110 
.02136 
.02306 


-1.8210 
-3.7309 
-5.5677 
-7.3981 


0.998 
.999 
.997 
.999 



important in order to correlate other solute properties for 
estimation purposes such as bioconcentration and soil 
adsorption. 

The advantages in using a log yf instead of log K o/w in 
defining H are: (1) solubility data are more readily 
available in the literature than K o/vf data, and (2) y£ is 
independent of any solvent/ water system. 

The utility of the gas chromatographic technique for 
measuring activity coefficients in organic solutions and 
the validity of the data obtained are now well estab- 
lished. The technique is particularly applicable to 
measurements at infinite dilution, the condition of 
primary interest and least accessible otherwise. There 
are several important advantages to this method for 
measuring y?: (D the speed with which the 
measurements are made, (2) the accuracy of the 
measurements, (3) the measurements are made at in- 
finite dilution, and <4) several solutes may be injected in- 
to the gas chromatograph simultaneously thus increasing 
the productivity rate. Although absolute y? value may be 
obtained by this method, the method is best suited for 
measuring values relative to some carefully studied stan- 
dards because of the relatively high volatility of octanol. 

In the octanol/water system at equilibrium the water 
phase is saturated with octanol(wo), and the octanol 
phase is saturated with water(ow). In order to derive an 
expression relating K 0/ r w to the solute aqueous solubility, 



we assume that y^*2s y 9 and y^° £ y^ where y£and y°. are 
the solute activity coefficients m pure water and pure oc- 
tanol, respectively. The extent to which these assump- 
tions are valid may be judged by the very good agree- 
ment between K / w values calculated via eq (10) and ex- 
perimental values (table 1) measured by a generator col- 
umn method [12,13]. 

The above method for determining K / w is best suited 
for volatile compounds. There are drawbacks to the 
method for relatively nonvolatile solutes. These com- 
pounds depending on their yg value could require a 
relatively long time to elute through the GC column, 
thus making the method impractical. 

For the cases where vapour pressure data are not 
available, Purnell has shown that carbon number (n c ) or 
boiling point {T b } may be substituted for log P°, i.e. 



hgK o/vi = k 1 n c + k 2 
where the k's are constants. 



(18) 
(19) 



Equation (19) is particularly useful when the 
homologous series cannot be described by the carbon 
numbers. The results listed in the table 2 show that there 
is an excellent correlation with these quantities. 

In summary it is clear that the chromatographic 
technique, particularly in association with the 
developments reported here, offers a primary route to 
rapid collection of large volumes of reliable values of 
J£ / w and other thermodynamic functions associated 
with the partitioning of organic material in the environ- 
ment. The infinite dilution method proposed in this 
paper has all the advantages of a chromatographic 
method in that there are no stringent demands placed on 
the purity of the solute or the amount of material re- 
quired to determine K 0/ - w . 
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